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Abstract This paper discusses the Lyapunov exponent A for small particles in a spatially and tempo- 
rally smooth flow in one dimension. Using a plausible model for the statistics of the velocity gradient in 
the vicinity of a particle, the Lyapunov exponent is obtained as a series expansion in the Stokes num- 
ber, St, which is a dimensionless measure of the importance of inertial effects. The approach described 
here can be extended to calculations of the Lyapunov exponents and of the correlation dimension for 
inertial particles in higher dimensions. It is also shown that there is correction to this theory which 
arises because the particles do not sample the velocity field ergodically. Using this non-ergodic correc- 
tion, it is found that (contrary to expectations) the first order term in the expansion of A does not 
vanish. 
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1 Introduction 

It is known that small particles suspended in complex flows, such as fully-developed turbulence, can 
cluster together p], so that they accumulate on a fractal attractor [2j. In cases where the flow has a 
compressible component, it is also possible for the trajectories of particles to coalesce [3J. These effects 
have been analysed by calculating the Lyapunov exponents of the particle trajectories (the role of 
Lyapunov exponents in characterising dynamical systems is discussed in [1]). The leading Lyapunov 
exponent A is the average of the logarithmic rate of divergence of particles with an infinitesimal 
separation 6x(t): 

~5x(t)~ 



A = lim — ( In 



Sx(0) 



(1) 



where (X) denotes the expectation value of X. If the leading Lyapunov exponent is negative, the 
particles aggregate together (the path coalescence effect) [5]. In cases where there is clustering onto a 
fractal attractor rather than aggregation, the dimension of the attractor can be estimated [B] from the 
full set of Lyapunov exponents of the particle trajectories by means of the Lyapunov dimension, given 
by the Kaplan- Yorke formula [7]. The clustering effect may play an important role in explaining the 
growth of rain drops [5] and of planetesimals [5], and deserves to be thoroughly understood. 
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Substantial progress has been made in analysing these clustering effects for a model with a random 
velocity field, in the limit where the correlation time of the velocity field is vanishingly small. All 
three Lyapunov exponents for the suspended particles were obtained as power series in a dimensionless 
parameter |10j . The methods used in that work involve a mapping to a Fokker- Planck equation, and 
appear to be restricted to the case of a vanishing correlation time, where the system is described by a 
Markovian stochastic process. This paper shows how the analysis developed in [10] may be extended 
to models where the velocity field has a finite correlation time, using a result from a companion paper, 
[TT] , This allows the theoretical analysis of clustering to be extended to physically realistic model 
velocity fields. 

It will be shown that there is an unexpected subtlety in the application of the finite correlation time 
model, which may complicate its application to models of turbulent velocity fields. For this reason, the 
problem is explored here in its simplest form, involving a one-dimensional flow, so that the difficulties 
can be explored in the simplest context. 

In order to analyse the clustering effect it is necessary to consider the dimensionless parameters of 
this problem. The Stokes number, a dimensionless measure of the inertia of the particles, is defined by 

St = — (2) 

7T 

where r is the correlation time for the fluctuations of the velocity gradient of the fluid, and 7 is the 
rate at which particle motion relative to the fluid is damped. The statistics of the velocity field are 
described by a dimensionless parameter known as the Kubo number: 

Ku=f (3) 

where u Q is a typical scale size of the velocity field and £ is a correlation length. In [TD] it is argued that 
for fully-developed turbulence the relevant length, time and velocity scales are those corresponding 
to the Kolmogorov microscales of the velocity field, and that dimensional considerations imply that 
the Kubo number is of order unity for turbulent flows. However, the Kubo number will remain as a 
parameter in any random flow model for turbulence. It is desirable to understand how the Lyapunov 
exponent varies as a function of both the Stokes number and the Kubo, so that the analytical properties 
of the approximations can be thoroughly understood. In section [3J the variables St and Ku will be 
replaced by variables to and k, which are dimcnsionally equivalent but which arc precisely defined. 

The model with vanishing correlation time which is analysed in [10] corresponds to considering the 
limit as Ku — ► 0. The application of models with with zero correlation time has been criticised [12 , 
because it is argued that that such models fail to get the leading order terms of series expansions for 
the Lyapunov exponents correct. In particular, it has been argued [12j that the difference between 
the fractal dimension D of the measure upon which particles cluster and the dimension d of the 
coordinate space is of second order in the Stokes parameter, that is d — D = 0(St 2 ), whereas for the 
zero correlation time model it is found that D — d = O(St) [TD]. Numerical evidence for the prediction 
that d — D = 0(St 2 ) is presented by Bee [13], who credits Balkovsky et al [14] with the theoretical basis 
for this result. This discrepancy is addressed by the model which is developed here, and it emerges that 
these issues can be quite subtle. For the one-dimensional model with zero correlation time it is known 
[5] that A/7 = O(St), and here it is shown that a simple application of the model with finite correlation 
time indicates that the leading order dependence of A/7 is 0(St 2 ), in line with the expectation that the 
finite correlation time model should have a different leading-order behaviour. A more careful analysis, 
however, leads to a surprising conclusion. The model with finite correlation time has the same leading 
order behaviour as the zero correlation time model, namely A/7 = O(St). Moreover, the prefactor is 
the same. 

Section [2] explains how the calculation of Lyapunov exponents is related to the problem discussed 
in the companion paper [TT] , which analyses the statistics of a class of stochastic differential equations 
which are driven by noise signal which has a finite correlation time. Section [3J shows how the series 
expansion for the expectation value derived in [11] is used to produce an estimate for the Lyapuonov 
exponent, which satisfies A/7 = (3(St 2 ). Section [4] shows that this simple estimate is based upon an 
unjustified assumption, namely that the positions of the particles sample the velocity gradient field 
ergodically. The leading order non-ergodic correction is computed, and it is shown that the correct 
leading order behaviour is A/7 = O(St). A different example of the importance of non-ergodic effects 
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in the dynamics of small particles in a random velocity field was commented upon in |15j . Section [5] is 
a summary of the conclusions, and points the way to applications in the systematic analysis of particle 
clustering in realistic two- and three-dimensional flows. 



2 A model for the Lyapunov exponent of particles in a random flow 

This section discusses the connection between the Lyapunov exponent and the equation of motion |T]) 
of the companion paper, [11J . Equations ((4]) to (|9|) follow arguments in [5] closely, but this material is 
included for the convenience of the reader. 

Consider a one-dimensional model for particles embedded in a fluid with velocity u(x,t). It is 
assumed that the motion of the particle relative to the fluid experiences a viscous damping force which 
is proportional to the velocity of the particle relative to the fluid. The equation of motion for a particle 
with position x and velocity v is therefore 

x = v , v = —"f[v — u(x, t)] . (4) 

The Lyapunov exponent characterises the time dependence of the infinitesimal separation Sx of a 
trajectory from a reference trajectory. To analyse this quantity, linearise the equation of motion to 
obtain 

Sv = Sx , Sv = — j[Sv — Y(t)Sx] (5) 

where 

Y(t) = ^(x(t),t) (6) 
is the velocity gradient at the current position of the particle. Now define a variable X: 

This has a simple relation to the Lyapunov exponent: 

A = (X) (8) 

and the equation of motion for X is 

X = -jX - X 2 + 7 Y(i) . (9) 

Thus the task of calculating A resolves into calculating the expectation value of a variable X which 
satisfies equation 

In turbulent or other unsteady, complex flows the velocity gradient Y(t) = du{x(t),t) / dx can be 
modelled as a stochastic variable. With this interpretation, equation © is equivalent to the stochastic 
differential equation ([I]) of the companion paper [11) . The velocity gradient has a finite correlation 
time, with a correlation function which depends upon the details of the flow. If the fluid velocity field 
u(x, t) is statistically homogeneous in space and time, it seems reasonable to model the fluctuations of 
the velocity gradient by a random variable which has mean value equal to zero. The companion paper 
discussed the case where the random variable Y(t) has a finite correlation time, by considering the 
following model for Y(t): 

Y = --Y + r)(t) (10) 

T 

with 

( v (t)) = o , (vihHh)) = 2D5{h - 1 2 ) . (ii) 

This is an Ornstein-Uhlenbeck process [16] , discussed in [17] . The correlation function of this process 
is 

{Y{h)Y{t 2 )) - DrexpHh - t 2 \/r) (12) 

which implies that 

° = ^r dt {Ymm (i3) 

so that the diffusion coefficient in (fTT]) is related to the correlation function of Y(t) by ([T5)l . 
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3 A series expansion for the Lyapunov exponent 

In the following equations ([9|), (fT0|) and (fTTj) will be transformed into a dimensionless form which is 
precisely the same as equations (4), (5) of the companion paper, The results developed there will 
be used to obtain a series expansion of the Lyapunov exponent, A = (X). It is convenient to transform 
to a new dimensionless time variable, t' = jt. And the variables X, Y are related to dimensionless 
variables by writing X = aX', Y = aY 1 , where the quantity a has dimensions of inverse time. By 
writing a = t^Dj the equations of motion , (|10[) can be written in the form 

¥L=Y'-X'- eX' 2 
dt' 

^ = _a;y' +wC (t') (14) 

where is a white noise with statistics 

(C(f)>=0, (C(ii)C(4))=25(ti-4). (15) 

Equations (|14|) . (|15|) are precisely the system which is analysed in the companion paper, [llj . The 
parameters e and uj in (|14p . and an additional parameter re, are defined by 

uj = — , k= Vlh^ . (16) 

7T 

The dimensionless parameters re and uj are precisely defined realisations of the dimensionless parameters 
introduced in section [TJ that is 

St = cj, Ku = re, e = KuSt 1/2 . (17) 

The Lyapunov exponent is 

A = a(X') = ^e(X') . (18) 

The companion paper gives a series expansion for the expectation value (X 1 ). Using (117|) to express e 
in terms of re and uj and using equation (30) of [llj gives 




2, ,2 



A = 7« w 



' +d(uj)K 2 UJ 2 + C2(uj)K i UJ i + 



(19) 



1 +CJ 

where the CAlS) are ratios of polynomials with the same degree, with C2(0) = 15 , 6*3(0) = 630 . . .. The 



leading order behaviour of the zero correlation time model is A ~ 7c 2 = jk 2 uj [5] . Thus equation (|19p 
is in agreement with the expectation that for a velocity field with finite correlation time, the leading 
order behaviour is different. It also has the satisfying feature that it is based upon a model which 
yields all of the coefficients of the series expansion exactly by using the operator algebra discussed in 
the companion paper. The technique appears to have a straightforward extension to multi-dimensional 
cases, and can also be adapted to provide the full series expansion of the subdominant Lyapunov 
exponents 10 and the correlation dimension [18] . 

This is however a deficiency in the calculation, which will be discussed in the next section. 



4 Non-ergodic correction to the velocity statistics 

Note that the model introduced in section [2] uses the seemingly uncontroversial assumption that the 
mean value of the velocity gradient Y(t) is zero. For a velocity field u(x,i) which is statistically 
homogeneous as a function of x, the expectation value of du/dx must certainly vanish. It is, however, 
possible for (Y(t)) to be non-zero because Y(t) is the velocity gradient evaluated at the current position 
of the particle, x(t). Because x(t) is itself dependent upon the history of the velocity field, Y(t) may 
not represent an ergodic sampling of this field. 

This non-ergodic correction may, in most circumstances, represent a small refinement of the theory. 
In the overdamped limit (as uj — > 0) it could, however, be of critical importance in determining the 
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leading order behaviour of A. In the limit as lu — > 0, the Ornstein-Uhlenbeck noise term Y' in equation 
(fT4f varies very slowly. The stochastic variable X 1 must therefore also vary very slowly, and the time- 
derivative in (|14p can therefore be neglected. In the overdamped limit, the equation of motion for X 1 
can therefore be approximated by 

Y' — X' — eX' 2 = . (20) 

This is a quadratic equation determining X' as a function of Y 1 . In the overdamped limit (where 
lu — > with k fixed), e = k^/lu — > 0. The solution of this quadratic equation is approximated by 
X' = Y' — eY' 2 + 0(e 2 ). Taking expectation values, and using (Ti"H|) , this implies that the Lyapunov 
exponent is 

A = 7e [(yv^' 2 >]+0( e 3 ). (2i) 

For the Ornstein-Uhlenbeck noise process in (fT4| . (Y 1 ) = and (Y' 2 } — lu, so that the approximate 
equation (|2ip gives A = je 2 ^ — jk 2 lu 2 , a result in precise agreement with the leading order term in 
(|18[) . If, however, a more refined approximation shows that (Y(t)) is non-zero, it may be discovered 
that the dominant contribution to A/7 does, in fact, come from a term which is linear in the Stokes 
number lu. 

In principle it is possible to expand (Y) as a series in lu and K. Only the leading term of this 
expansion is considered here. In the limit where lu <C 1 a particle is simply advected by the flow and 
the equation of motion may be approximated by the advection equation, x — u(x(t),t). The position 
of the particle is then 

(22) 



(23) 



x(t)=x(0) + dt' u(x(t'),t') + 0{lu) . 
Jo 

Now use this approximation to calculate the time-average of Y(t), that is 



(Y(t)) = lim 

1 — >oc 







T , du , , . . 
dt fo(x(t),t) 



Using equation (f2"2")l . 







dt 



du 
dx 



(x(0),t) 



d 2 u { 
dx 2 



x(0),t) / dt' u(x(0),t') 



(24) 



Averaging this, using translational invariance and assuming T>t gives 



(Y(t))= lim I / T dt Tdt' (^tt(0 s t)«(a!, «'))+.. 
dt (flgu(0,t)«(0,Q)) + ... . 



(25) 



The correlation function in (I25|) is easily related to that in (|13| . To this end define correlation functions 
C (t), Ci(t) and C 2 (t) as follows: 

}2„, 



C (t) 



_9_ 

<9.i: 



du 
dx 



(x, t)u{x, 0) 



Ci(t) 



(i,t)«(i,0) , C 3 (*) = 



9u . . 9?i . 



(26) 



Note that Co(t) = C'i(t) + C*2(t). Consider the evaluation of Co(t) by averaging over the spatial 
coordinate: because the quantity being averaged is a spatial derivative, 



C (t) 



lim — 

L — >oo L 



dx 



0_ 

dx 



du 
dx 



(x, t)u(x, 0) 



= 



(27) 



implying that Ci(t) 
tion of D in (fT3|). 



-Czit). Now return to the evaluation of (Y(t)) using ([25|) . Recalling the defini- 



dt C*i(t) 



dt C 2 (t) 



-t D 



(28) 



From this it follows that (Y') = Ty/D/^j = e, and hence, using (f2"Tj) . one obtains the following estimate 
for the Lyapunov exponent 

A = - 7 e 2 = -jk 2 lu . (29) 
This is exactly the form obtained using the zero correlation time model, discussed in [5]. 
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5 Discussion 

This paper has discussed a new approach to extending the theoretical understanding of the dynamics 
of particles suspended in random flows, through making series expansions of the Lyapunov exponent 
in powers of the Stokes number. Until now, this has only been possible for models of the fluid velocity 
field with zero correlation time. This paper has shown how the techniques for obtaining the series 
expansion can be extended to a class of models for flows with a finite correlation time. 

It has been questioned whether the leading order behaviour of these series is correct for the zero 
correlation time model where, k — > 0. For example, it has been proposed that the difference between 
the fractal dimension of the attractor and the space dimension is 0(St 2 ) for real flows, whereas it is 
found to be O(St) for the zero correlation time model. A direct application of the finite correlation 
time model, as considered in section |4l appears to be consistent with this expectation, in that the 
leading order dependence of A/7 appears to be 0(St 2 ) for the finite correlation time case, whereas it 
is O(St) for the zero correlation time model. 

However, a more careful examination of the problem leads to a different conclusion. The subtlety 
lies in considering the statistics of the velocity gradient of the particle at its current position, Y(t) = 
du/dx(x(t),t). In the case where the correlation time approaches zero, statistics of Y(t) are clearly 
that same as those of du/dx sampled at a randomly chosen point in space. When the correlation 
time of the velocity field is finite, this ergodic relation can be questioned, because the position x(t) 
has been influenced by the history of the velocity field, which is correlated with its current value. In 
section [4] it was shown that these correlations are significant, and imply that the mean value of the 
velocity gradient at the particle position Y(t) is non-zero, despite the fact that the mean value of 
du/dx is obviously zero. When this correction is accounted for, the leading order behaviour of the 
finite correlation time model becomes exactly the same as the that of the zero correlation time model, 
including the prefactor. 

The principal motivation for this study was to develop a framework for extending the analysis of 
the dynamics of small particles in a random flow to models with finite correlation time. The approach 
developed in the companion paper and in section [3] of this work is well suited to this end. The result in 
section [4] indicate that the statistics of the velocity gradient field as sampled by the particles have to be 
modelled carefully in order to obtain precise results. These conclusions do not contradict the numerical 
results on the fractal dimension presented in [13j , because the calculation of the fractal dimension is a 
different problem from that of the largest Lyapunov dimension, considered here. But they do indicate 
that the extension to modelling physically realistic flows in two and three dimensions will be a delicate 
issue. 

Acknowledgement. I have benefitted from extensive discussions with Bernhard Mehlig, and the 
calculation in section |4] was prompted by his numerical experiments giving evidence for the result 
reported there. I am also grateful to Michael Morgan for helpful comments on this paper. 
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